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I. INTRODUCTION 

One of the most attractive aspects of the vacuum of 
Quantum Chromodynamics (QCD), the theory of strong 
interactions, is its non-pcrturbative nature. Mainly by 
means of Lattice QCD simulations at zero quark chem- 
ical potential [1-5] it is established that two crossovers 
take place in a narrow range of temperature; one for 
quark deconnnement, and another one for the (approx- 
imate) restoration of chiral symmetry. Besides, power- 
ful analytic and semi-analytic techniques have been de- 
veloped to understand the coupling between chiral sym- 
metry restoration and deconfmement, see 0, 0] and ref- 
erences therein. Moreover, Lattice QCD and Operator 
Product Expansion (OPE in the following) of the corre- 
lators of hadronic currents show that the QCD vacuum 
can be characterized by several quark, gluon and mixed 
condensates @. 

A fruitful theoretical approach to the physics of strong 
interactions, which is capable to capture some of the non- 
pcrturbative properties of the QCD vacuum, is the use 
of effective chiral models. One of them is the Nambu- 
Jona-Lasinio (NJL) model @ (see Refs. for reviews), 
in which the QCD gluon-mediated interactions are re- 
placed by effective interactions among quarks, which are 
built in order to respect the global symmetries of QCD. 
Under some approximations, it is possible to derive the 
NJL model effective interaction kernel from first princi- 
ples QCD, see [ll|, E3- Besides, a linear a— model aug- 
mented with a Yukawa- type coupling to quarks, named 
Quark-Meson model (QM model in the following) , has 
been developed as an effective model of QCD ID, [lj] . In 
this model, quartic meson self-couplings allow to absorb 
the cutoff dependence of the coupling constants; more- 
over, following an idea by Weinberg 15|. tree level prop- 



agating mesons are sufficient to avoid triviality of the 
NJL model in 3 + 1 dimensions in the one-loop approx- 
imation [T|| (see [lH for an excellent discussion about 
these points). 

The chiral models are widely used to map qualitatively, 
and to some extent also quantitatively, the phase diagram 
of strongly interacting matter along several directions like 
temperature, chemical potential, isospin chemical poten- 
tial and external fields [j"7r!30l]. It is thus interesting to 
use these models to compute other quantities, which are 
related to the QCD vacuum condensates. As a matter of 
fact, the chiral models allow for a self-consistent treat- 
ment of the spontaneous chiral symmetry breaking in 
the vacuum. Therefore, after the ground state properties 
(i.e., the chiral condensate in the NJL model, or the ex- 
pectation value of the cr-field in the quark-meson model) 
under external factors are computed, it is straightforward 
to estimate numerical values of other condensates, using 
the one-loop propagators of the theory. 

It has been realized that external fields can induce 
QCD condensates that are absent otherwise [3l|. Of 
particular interest for this article is magnetic moment, 
(fE^f) where / denotes the fermion field of the flavor 



f-th, and W 
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one can write, according to [311 ] . 

(fWf)= X (ff)Qf\eB\ 



(1) 
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and x is a constant independent on flavor, which is 
dubbed magnetic susceptibility of the quark condensate. 
In [3l[ it is proved that the role of the condensate (HJ 
to QCD sum rules in external fields is significant, and it 
cannot be ignored. The quantity y has been computed 
by means of special sum rules [3ll435l |. OPE combined 
with Pion Dominance (3(| , holography [13, [38| , instan- 
ton vacuum model [39j . analytically from the zero mode 
of the Dirac operator in the background of a SU (2) in- 
stanton [40| . and on the Lattice in two color quenched 
simulations at zero and finite temperature |4lj . It has 
also been suggested that in the photoproduction of lep- 
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ton pairs, the interference of the Drell-Yan amplitude 
with the amplitude of a process where the photon couples 
to quarks through its chiral-odd distribution amplitude, 
which is normalized to the magnetic susceptibility of the 
QCD vacuum, is possible [42|. This interference allows 
in principle to access the chiral odd transversity parton 
distribution in the proton. Therefore, this quantity is in- 
teresting both theoretically and phcnomenologically. The 
several estimates, that we briefly review in Section III, 
lead to the numerical value of x as follows: 



X(ff) = 40 - 70 MeV 



(2) 



A second quantity, which embeds non-linear effects at 
large fields, is the polarization, fit, defined as 



(//> 



E f = </£ 12 /) , 



(3) 



which has been computed on the Lattice in [41| for a 
wide range of magnetic fields, in the framework of two- 
color QCD with quenched fermions. At small fields 
A 4 / — \xQf £ B\ naturally; at large fields, non-linear ef- 
fects dominate and an interesting saturation of pf to the 
asymptotic value p^ = 1 is measured. According to [4l[ 
the behavior of the polarization as a function of eB in 
the whole range examined, can be described by a sim- 
ple inverse tangent function. Besides, magnetization of 
the QCD vacuum has been computed in the strong field 
limit in (43j using perturbative QCD, where it is found 
it grows as BlogB. 

In this article, we compute the magnetic susceptibil- 
ity of the quark condensate by means of the NJL and 
the QM models. This study is interesting because in the 
chiral models, it is possible to compute self-consistently 
the numerical values of the condensates as a function 
of eB, once the parameters are fixed to reproduce some 
characteristic of the QCD vacuum. We firstly perform a 
numerical study of the problem, which is then comple- 
mented by some analytic estimate of the same quantity 
within the renormalized QM model. Moreover, we com- 
pute the polarization of quarks at small as well as large 
fields, both numerically and analytically. In agreement 
with the Lattice results [4l|, we also measure a satura- 
tion of /i/ to one at large fields, in the case of the effective 
models. Our results push towards the interpretation of 
the saturation as a non-artifact of the Lattice. On the 
contrary, we can offer a simple physical understanding of 
this behavior, in terms of lowest Landau level dominance 
of the chiral condensate. As a matter of fact, using the 
simple equations of the models for the chiral condensate 
and for the magnetic moment, we can show that at large 
magnetic field pf has to saturate to one, because in this 
limit the higher Landau levels are expelled from the chi- 
ral condensate; as a consequence, the ratio of the two 
approaches one asymptotically. 

We also obtain a saturation of the polarization within 
the renormalized QM model. There are some differ- 
ences, however, in comparison with the results of the non- 
renormalizcd models. In the former case, the asymptotic 



value of pf is charge-dependent; moreover, the interpre- 
tation of the saturation as a lowest Landau level (LLL) 
dominance is not straightforward, because the renormal- 
ized contribution of the higher Landau levels is important 
in the chiral condensate, even in the limit of very strong 
fields. It is possible that the results obtained within the 
renormalized model are a little bit far from true QCD. 
As a matter of fact, in the renormalized model we as- 
sume that the quark self-energy is independent on mo- 
mentum; thus, when we take the limit of infinite quark 
momentum in the gap equation, and absorb the ultravio- 
let divergences by means of counterterms and renormal- 
ization conditions, we implicitly assume that that quark 
mass at large momenta is equal to its value at zero mo- 
mentum. We know that this is not true, see for exam- 
ple [iil . |45| : even in the renormalized theory, the quark 
self-energy naturally cuts off the large momenta, lead- 
ing to LLL dominance in the traces of quark propagator 
which are relevant for our study. Nevertheless, it is worth 
to study this problem within the renormalized QM model 
in its simplest version, because it helps to understand the 
structure of this theory under the influence of a strong 
magnetic field. 



In our calculations we neglect, for simplicity, the pos- 
sible condensation of p— mesons at strong fields [H, 
HtJ- Vector meson dominance [46| and Sakai-Sugimoto 
model [48| suggest for the condensation a critical value 
of eB c w m 2 sa 0.57 GeV 2 , where m p is the p— meson 
mass in the vacuum. Beside these, a NJL-based calcu- 
lation within the lowest Landau level (LLL) approxima- 
tion [47} predicts p— meson condensation at strong fields 
as well, even if in the latter case it is hard to estimate 
exactly eB Cl mainly because of the uncertainty of the 
parameters of the model. It would certainly be interest- 
ing to address this problem within our calculations, in 
which not only the LLL but also the higher Landau lev- 
els are considered, and in which the spontaneous breaking 
of chiral symmetry is kept into account self-consistently. 
However, this would complicate significantly the calcula- 
tional setup. Therefore, for simplicity we leave this issue 
to a future project. 



The plan of the article is as follows. In Section II we 
describe the QM and NJL models and fix our notation. 
In Section III we discuss our numerical results for the 
polarization of quarks, and compute the magnetic sus- 
ceptibility of the quark condensate. In Section IV we 
compute the renormalized Quantum Effective Potential 
(QEP) of the QM model in a magnetic background, in 
the one-loop approximation, and compute analytically 
the solution of the gap equation in the weak field case, 
and semi-analytically in the strong field limit. We then 
use the results to estimate pf and x- Finally, in Section 
V we draw our conclusions. 
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II. CHIRAL MODELS COUPLED TO A 
MAGNETIC FIELD 

In the first part of this article, we derive numerical 
results for the spin polarization of quarks in a magnetic 
field, and for the magnetic susceptibility of the quark con- 
densate, using two effective chiral models: the Nambu- 
Jona-Lasinio (NJL) model, and the Quark-Meson (QM) 
model. Wc describe the two models in some detail in 
this Section. We work in the Landau gauge and take the 
magnetic field along the z— axis, B = (0, 0, B). 

A. Quark-Meson model 

In the QM model, a meson sector described by the lin- 
ear sigma model lagrangian, is coupled to quarks via a 
Yukawa-type interaction. The model is renormalizable in 
D = 3 + 1 dimensions. However, since we adopt the point 
of view of it as an effective description of QCD, it is not 
necessary to use the renormalizcd version of the model 
itself. On the contrary, it is enough to fix an ultravio- 
let scale to cutoff the divergent expectation values; the 
UV scale is then chosen phenomenologically, by requiring 
that the numerical value of the chiral condensate in the 
vacuum obtained within the model, is consistent with the 
results obtained from the sum rules (49j . This is a rough 
approximation of the QCD effective quark mass, which 
smoothly decays at large momenta [Ij, [45[ . In Section 
IV we will use a renormalized version of the model, to 
derive semi-analytically some results in the two regimes 
of weak and strong fields. 

The lagrangian density of the model is given by 

C = q [iD^ 1 ~ 9{a + il5i~ ■ tt)] q 

+ \{d,a) 2 + l -(d^) 2 -U(a^) . (4) 

In the above equation, q corresponds to a quark field 
in the fundamental representation of color group SU(S) 
and flavor group SU(2); the covariant derivative, = 
9fi — QfeAfj,, describes the coupling to the background 
magnetic field, where Qf denotes the charge of the flavor 
/. Besides, a, 7r correspond to the scalar singlet and the 
pseudo-scalar iso-triplet fields, respectively. The poten- 
tial U describes tree-level interactions among the meson 
fields. In this article, we take its analytic form as 

U(a,Tv) = ^{a 2 + iv 2 -v 2 f -ha , (5) 

where the first addendum is chiral invariant; the second 
one describes a soft explicit breaking of chiral symmetry, 
and it is thus responsible for the non-zero value of the 
pion mass. For h = 0, the interaction terms of the model 
are invariant under SU{2) V ® SU{2) A ® U(l) v . This 
group is broken explicitly to U(l) v ® U{1)\ ® E/(l)y if 
the magnetic field is coupled to the quarks, because of 
the different electric charge of u and d quarks. Here, the 



superscript 3 in the V and A groups denotes the trans- 
formations generated by T3, T375 respectively. There- 
fore, the chiral group in presence of a magnetic field is 
U(l)y (8 This group is then explicitly broken by 

h-tevm to U(l) v . 

In this article, we restrict ourselves to the large-A^ 
(that is, one-loop) approximation, which amounts to con- 
sider mesons as classical fields, and integrate only over 
fcrmions in the generating functional of the theory to ob- 
tain the Quantum Effective Potential (QEP). As a matter 
of fact, quantum corrections arising from meson bubbles 
are suppressed of a factor 1/N C with respect to case of 
the fermion bubble. In the integration process, the me- 
son fields are fixed to their classical expectation value, 
(tv) = and (a) ^ (in particular, a has the quantum 
numbers of the chiral condensate, (qq)). The physical 
value of (a) will be then determined by minimization of 
the QEP. 

To compute QEP in presence of a magnetic back- 
ground, we use the Leung-Ritus-Wang method [5(| which 
allows to write down the quark propagator for the flavor 
/ in terms of Landau levels, 

*/(*.»> = tj dJ ^r lE ^p:^M^ , 

(6) 

where Ep(x) corresponds to the eigenfunction of a 
charged fermion in magnetic field, and Ep(x) = 
lo{Ep(x)) T70. In the above equation, 

P={po,0,Qyj2k\Q f eB\,p 3 ) , (7) 

where k = 0,1,2,... labels the k th Landau level, and 
Q = sign(Qf), with Qf denoting the charge of the fla- 
vor /; Afe is a projector in Dirac space which keeps into 
account the degeneracy of the Landau levels; it is given 
by 

A fc = S k0 [V+S Q , +1 + V-6q,-i] + (1 - 6 k0 )I , (8) 

where V± arc spin projectors and I is the identity matrix 
in Dirac spinor indices. At the one-loop level, the QEP 
then reads 

V=^{<j 2 -v 2 ) 2 -ha 

/ k J 

In the above equation we have defined 

uj k (p 3 ) = ^ P l + 2k\Q f eB\+ml , (10) 

with TTtq — ga; (3 k = 2 — 5 k o counts the degeneracy of the 
Landau levels. 

The one-loop fermion contribution, which corresponds 
to the last addendum in the r.h.s. of Eq. ([9]), is divergent 
in the ultraviolet. In order to regularize it, we adopt a 
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smooth regulator Ua as in [26| . which is more suitable, 
from the numerical point of view, in our model calcu- 
lation with respect to the hard-cutoff which is used in 
analogous calculations without magnetic field. In this 
article we chose 



U A 



A 



2iV 



A 2W + (p2 + 2\Q f eB\k) 



N 



N = 5 



(11) 



The (more usual) 3-momentum cutoff regularization 
scheme is recovered in the limit N — > oo; we notice that, 
even if the choice N = 5 may seem arbitrary to some 
extent, it is not more arbitrary than the choice of the 
hard cutoff scheme, that is, of a regularization scheme. 
In effective models, the choice of a regularization scheme 
is a part of the definition of the model itself. Momentum 
integrals are understood as follows: 

2>/|-W£*- ,12) 

n J n J 

Once the expectation value of a is computed as a func- 
tion of eB by a procedure of minimization of V, we com- 
pute the expectation values that are relevant to the con- 
text. To begin with, we consider the chiral condensate 
for the flavor /, (//) = —Tr[Sf(x,x)], with Sf given by 
Eq. ([5]) . It is straightforward to derive the relation 



(//) 



-N, 



\QfeB\ 
2ir 



oo . 

w 

k=0 J 



dp3 



2tt uj k (P3) 



(13) 



where the divergent integral on the r.h.s. of the above 
equation has to be understood regularized as in (|T2"1) . 
From Eq. (fT5|) we notice that the prescription (|I2[) is al- 
most equivalent to the introduction of a running effective 
quark mass, 



= gaQ (A 2 - p 2 3 - 2k\Q f eB\) 



(14) 



that can be considered as a rough approximation to the 
effective running quark mass in QCD (45j which decays 
at large quark momenta, see also the discussion in (4J). 
Once the scale A is fixed, the Landau levels with n > 1 
are removed from the chiral condensate if eB 3> A 2 . 
Next we turn to the magnetic moment for the flavor /, 



where 



(fWf) =-Tr[W S f (x,x)] 



(15) 



(16) 



is the relativistic spin operator. We take B = (0,0, B); 
in this case, only E 12 = £/ is non- vanishing. Using the 
properties of 7— matrices it is easy to show that only 
the Lowest Landau Level (LLL) gives a non-vanishing 
contribution to the trace: 



Nr. 



Qf\eB\ f dp 3 m q 



2tt 



2tt uj Q (p 3 ) 



(17) 



where ujq = uj^o- Once again, the divergent integral on 
the r.h.s. of the above equation has to be understood 
regularized following the prescription in Eq. (JT2J - 

For the QM model, the parameters are tuned as fol- 
lows. We fix the constituent quark mass in the vacuum to 
a phenomenological value, m q = 335 MeV; furthermore, 
we require that in the vacuum the following condition 
holds: 



dV(a, B = 0) 



da 







(18) 



which implies that (a) = / w in the vacuum, thus m q = 
gfir- As a consequence we find g = 3.62. The parameter 
h is fixed by the condition h — /^m 2 , where is the 
average value of the pion mass in the vacuum and / w = 
92.4 MeV. We have then h = 0.047m?. To determine v 
and A we solve simultaneously Eq. ([T8j and 



d 2 V 



da 2 



(19) 



with m a = 700 MeV. The divergences in these two 
equations are cured with the prescription (|12l) . In the 
UV regulator wc chose A = 560 MeV which implies 
(uu) = (-231 MeV) 3 . This results in the values A = 4.67 
and v 2 = —1.8m 2 . Before going ahead, we notice that 
in the non-renormalized model, the 1-loop fermion con- 
tribution regularized at the scale A is included into the 
conditions JT8l) and (IT9l). 



B. Nambu-Jona-Lasinio model 

The quark lagrangian density of the NJL model is given 

by 

£ = q(i^D tl -m )q + G[(qq) 2 + {iq l5 Tq) 2 } ; (20) 

here q is the quark Dirac spinor in the fundamental repre- 
sentation of the flavor SU{2) and the color group; r cor- 
respond to the Pauli matrices in flavor space. A sum over 
color and flavor is understood. Once again, the covariant 
derivative embeds the QED coupling of the quarks with 
the external magnetic field. The explicit soft breaking of 
chiral symmetry in this model is achieved by introducing 
a current quark mass, mo, whose numerical value is fixed 
by fitting the vacuum pion mass. The NJL model is not 
renormalizable in the usual sense for D > 2 space-time 
dimensions; it is renormalizable in D < 4 in the one-loop 
approximation. At D = 4 it might represent a trivial 
theory of non-interacting bosons if the ultraviolet cut- 
off is allowed to be infinite [13, (3 ■ However, adopting 
the point of view of the NJL model as an effective (and 
rough) description of low-energy QCD, it is not necessary 
to require the cutoff to be infinite: the cutoff can be inter- 
preted as a physical quantity which makes the effective 
quark mass to be constant at small quark momenta, and 
vanishing at large momenta, thus mimicking roughly the 
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running effective mass of real QCD 45|. It is worth to 
notice that a non-local interaction among quarks with the 
same quantum numbers of the four-fermion local term in 
Eq. (p?0|) can be derived by the QCD action under some 
approximation, taking the low energy limit of the lat- 
ter 0, [H| . For simplicity, we treat here only the case of 
local interaction. 

Once again, the one-loop fermion contribution can be 
obtained within the Leung-Ritus-Wang method: 



V = GE 2 - N c 

/=«,< 



\QfeB\ 
2tt 



7— oo -* n 



(21) 



with ui n (p) - 
2GE with E 



p* + 2n\Q f eB\ + M 2 , and M q = m 

uu + dd) . Equations (fT5|) and (|15p are 
still valid in the NJL model, with the replacement m q — > 
Mq. A comparison of the QEPs of the QM model and 
of the NJL model shows that at the one-loop level, the 
two models differ only for the classic part of the effective 
meson potential, and for the definition of the constituent 
quark mass. 

The parameters in the NJL model are fixed as follows. 
The bare quark mass is computed by virtue of the Gell- 
Mann-Oakes-Renner relation, m 2 . f 2 = — 2mo(uu), which 
is satisfied in the NJL model [1(|- Then, two equations 
are solved simultaneously: one for f„ and one for the 
chiral condensate in the vacuum at zero magnetic field 
strength. Requiring that f n = 92.4 MeV and (uu) = 
(-253 MeV) 3 , this procedure gives the numerical values 
of G and A. Thus we find A = 626.76 MeV, G = 2.02/A 2 
and mo = 5 MeV. 



III. NUMERICAL RESULTS 

In this Section, we collect our numerical results which 
are relevant for the computation of the magnetic suscep- 
tibility of the chiral condensate, and of the spin polar- 
ization. From the numerical point of view, it is more 
convenient to compute firstly the latter; then, a fit of the 
polarization data at small fields will enable to extract the 
value of the magnetic susceptibility of the quark conden- 
sate. For what concerns data about chiral condensate, 
magnetic moment and polarization, we plot results only 
for the QM model, since the results obtained within the 
NJL model are qualitatively very similar to those ob- 
tained within the former model. Then, we will give the 
final result for the chiral magnetization for the two mod- 
els considered here. 



A. Polarization 

The physical value of the (total) chiral condensate E 
for the NJL model, and of a for the QM model, are ob- 
tained numerically by a minimization procedure of the 
QEP (f2"Tj) for any value of eB. Then, we make use of 
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FIG. 1. Upper panel. Chiral condensates of M-quarks (red) 
and d-quarks (blue), in units of the same quantities at zero 
magnetic field, as a function of the magnetic field. Lower 
panel. Expectation value of the magnetic moment operator, 
in units of f%, as a function of eB. Data correspond to the 
QM model. 



Eqs. (fT3"|) and (fT5"|) . with the replacement (fT2")l , to com- 
pute the chiral condensate and the magnetic moment for 
each flavor. 

In the upper panel of Fig. [TJ we plot the chiral con- 
densates for u and d quarks, as a function of eB, for the 
QM model. The magnetic field splits the two quantities 
because of the different charge for the two quarks. The 
small oscillations, which are more evident for the case of 
the it-quark, are an artifact of the regularization scheme, 
and disappear if smoother regulators are used, see the 
discussion in [2^. In the regime of weak fields, our data 
are consistent with the scaling (//) oc \eB\ 2 /M where 
M denotes some mass scale; in the strong field limit we 
find instead (//) oc |eB| 3 / 2 . The behavior of the quark 
condensate as a function of magnetic field is in agreement 
with the magnetic catalysis scenario [13, HH . 

In the lower panel of Fig. [T] we plot our data for the 
expectation value of the magnetic moment. Data corre- 
spond to the QM model (for the NJL model we obtain 
similar results). At weak fields, E/ oc \eB\ as expected 
from Eq. (|15p . In the strong field limit, non-linearity 
arise because of the scaling of quark mass (or chiral con- 
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FIG. 2. Polarization of it-quarks (red) and d— quarks (blue) 
as a function the magnetic field strength, for the QM model. 



densate); we find £/ oc |e£?| 3 / 2 in this limit. 

In Fig. [2] we plot our results for the polarization. 
Data are obtained by the previous ones, using the defi- 
nition ©. At small fields, the polarization clearly grows 
linearly with the magnetic field. This is a natural conse- 
quence of the linear behavior of the magnetic moment as 
a function of eB for small fields, see Fig. [JJ On the other 
hand, within the chiral models we measure a saturation of 
fif at large values of eB, to an asymptotic value poo = 1. 
This conclusion remains unchanged if wc consider the 
NJL model, and it is in agreement with the recent Lat- 
tice findings [4l|. It should be noticed that, at least for 
the u— quark, saturation is achieved before the expected 
threshold for p— meson condensation (4fj| - l48l |. Therefore, 
our expectation is that our result is stable also if vector 
meson condensation is considered. 

The saturation to the asymptotic value /ioo = 1 of po- 
larization is naturally understood within the models we 
investigate, as a LLL dominance in the chiral condensate 
(i.e., full polarization). As a matter of fact, and (//) 
turn to be proportional in the strong field limit, since 
only the LLL gives a contribution to the the latter, com- 
pare Eq. and (fTT)) which imply 



Hf = l- 



(//) 



HLL 



(//> 



(22) 



where (//)hll corresponds to the higher Landau lev- 
els contribution to the chiral condensate. In the strong 
field limit (//)hll — > because of Eq. (fT2"j) which is a 
rough approximation to the QCD effective quark mass, 
as depicted in Eq. hence, pf has to approach the 

asymptotic value fioo ~ 1. On the other hand, in the 
weak field limit (//)hll — > (//) & n d the proportionality 
among Y,f and (//) is lost. 

In Rcf. [4l|, data of polarization are fit by means of 
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FIG. 3. Upper panel. Polarization data for it— quark (grey 
dots) and fitting curves. See the text for details. Lower panel. 
Polarization data and fitting functions for d— quark. Results 
correspond to the QM model. 



three different functions, namely 



clas _ 

fif — MOO 



quan 

M f — Moo 



coth 



Mo 



2 coth 



3xQfeB_ 

Moo 

ZxQfeB X Q f eB 



3 X Q f eB 
- coth 



trig 

fj, f = arctan 



Moo 

irxQfeB 



2/i c 



(23) 

, (24) 
(25) 



The three functions above share the behavior at the ori- 
gin, n ,f m \xQ/eB\, and the asymptotic one, /i/ ~ Hoc- 
In [4l| a two-parameter fit on x an d Moo is performed; the 
model with the lowest chi-squared per degrees of freedom 
is represented by the trigonometric one, Eq. (|2"5")) . In- 
spired by these results we have tried to fit our data using 
the same functions ((23|) - ((25|) . In our model calculation, 
the asymptotic value fioo = 1 is achieved straightfor- 
wardly, therefore it is enough to perform a one-parameter 
fit leaving x as a free parameter. The results of this pro- 
cedure are collected in Fig. [3j where we plot our data 
of quark's polarization as a function of eB (gray dots) 
and the three fitting functions with /i^ = 1. From Fig. [3] 
we read that both the trigonometric and the classic func- 
tions do not adapt well to our data (for the NJL model we 
obtain similar results). In particular, both of the afore- 
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mentioned functions overestimate x, and reach slowly the 
asymptotic value. On the other hand, the quantum func- 
tion in Eq. ((24)) offers a better description of our data in 
the whole range of eB examined here. We have checked 
that the value of \ obtained within this fit overestimates 
the value obtained by a weak-field linear fit only of the 
10%, see below. Moreover, the fit functions smoothly fol- 
low the large eB data to the asymptotic value fXf = 1. 
We conclude that within our model, the quantum fitting 
function is a more faithful representation of data. The 
difference with Lattice data is probably due to the fact 
that the latter ones are obtained with quenched (hence, 
non-dynamical) fermions, while in our case dynamical 
fermions are considered. It will be interesting, therefore, 
to compare the model predictions with Lattice data ob- 
tained with other kinds of fermions in the near future. 

Before going ahead, it is interesting that our data on 
polarization, and our interpretation of the saturation of 
the latter at large fields, gives a quantitative estimate of 
the goodness of the LLL approximation for the models at 
hand. In particular, at eB c ps 0.6 GeV 2 , which in turns 
is the estimated critical value for vector meson conden- 
sation [46|, [48| , the LLL almost dominates the chiral con- 
densate of the u— quark; on the other hand, at this value 
of eB the LLL dominance has not yet been achieved for 
the d— quark condensate; but the higher Landau levels at 
this value of eB give a contribution almost to the 20% 
of the chiral condensate. This is an indirect check of the 
consistence of the LLL approximation used in |47j , which 
in turn should be good within a 20% accuracy. 

B. Magnetic susceptibility of the quark condensate 

At small fields fj,f = \xQ/eB\ from Eq. ([1]). Hence, 
we use the data on polarization at small fields, to obtain 
the numerical value of the magnetic susceptibility of the 
chiral condensate. Our results are as follows: 

X « -4.3 GcV~ 2 , NJL (26) 
X « -5.25 GcV~ 2 , QM (27) 

respectively for the NJL model and the QM model. To 
obtain the numerical values above we have used data for 
eB up to 5m 2 « 0.1 GeV 2 , which are then fit using 
a linear law. Using the numerical values of the chiral 
condensate in the two models, we obtain 

X(ff) « 69 McV , NJL (28) 
X{ff) « 65 MeV , QM (29) 

The numerical values of x that we obtain within the ef- 
fective models are in fair agreement with recent results, 
see Table I. In our model calculations, the role of the 
rcnormalization scale is played approximately by the ul- 
traviolet cutoff, that is A in Eq. (jlip . which is equal to 
0.560 GeV in the QM model, and 0.627 GeV in the NJL 
model. 

To facilitate the comparison with previous estimates, 
we review briefly the frameworks in which the results in 



Table I are obtained. In [36[ the following result is found, 
within OPE combined with Pion Dominance (we follow 
when possible the notation used in [ill)): 

X FD = — c x ' ' 2 , Pion Dominance (30) 

OTT r ^ 

with F n = VZfn = 130.7 McV and c x = 2; the estimate 
of [36j is done at a renormalization point M — 0.5 GeV. 
It is remarkable that Eq. (f3T)]) has been reproduced re- 
cently within AdS/Q CD approach in 38). Probably, this 
is the result more comparable to our estimate, because 
the reference scales in (3(| and in this article are very 
close. Within our model calculations we find c^ JL = 1.93 
and c^ M = 2.36. Using the numerical value of F„ and 
c x we get x PD = —4.45 GeV~ 2 , which agrees within 3% 
with our NJL model result, and within 18% with our QM 
model result. 

In [37| the authors find c x = 2.15 within hard- wall 
holographic approach, at the scale M <C 1 GeV. The 
results of 1371 are thus in very good parametric agree- 
ment with [36j; on the other hand, the numerical value 
of in the holographic model is smaller than the one 
used in [36| . pushing the holo grap hic prediction for x to 
slightly higher values than in [361 ) . However, the scale at 
which the result of [37l ) is valid should be much smaller 
than M = 1 GeV, thus some quantitative disagreement 
with [36| is expected. As the authors have explained, it 
might be possible to tunc the parameters of the holo- 
graphic model, mainly the chiral condensate, to repro- 
duce the correct value of F^: their numerical tests sug- 
gest that by changing the ratio (//)/m p of a factor of 
8, then the numerical value of c x is influenced only by 
a 5%. It is therefore plausible that a best tuning makes 
the quantitative prediction of [13| much closer to the es- 
timate of [36| . 

In [3^ ] an estimate of x within the instanton vacuum 
model has been performed beyond the chiral limit, both 
for light and for strange quarks (the result quoted in Ta- 
ble I corresponds to the light quarks; for the strange 
quark, Xs/Xu.d ~ 0.15 is found). Taking into account 
the numerical value of the chiral condensate in the in- 
stanton vacuum, the numerical estimate of [39[ leads to 
X = -2.5 ± 0.15 GeV" 2 at the scale M = 1 GeV. An 
analytic estimate within a similar framework has been 
obtained in [i(| , in which the zero- mode of the Dirac op- 
erator in the background of a SU(2) instanton is used 
to compute the relevant expectation values. The result 
of [H gives x = -3.52 GcV~ 2 at M « 1 GeV. 

In (41j the result x = —1-547 GcV~ 2 is achieved 
within a two-color simulation with quenched fermions. 
It is interesting that in [il| the same quantity has been 
computed also at finite temperature in the confinement 
phase, at T = 0.82T C , and the result seems to be in- 
dependent on temperature. The reference scale of [il| . 
determined by the inverse lattice spacing, is M sa 2 GeV. 
Therefore the lattice results are not quantitatively com- 
parable with our model calculation. However, they share 
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an important feature with the results presented here, 
namely the saturation of the polarization at large values 
of the magnetic field. Finally, estimates of the magnetic 
susceptibility of the chiral condensate by means of sev- 
eral QCD sum rules there exist (3ll - [35j . The results are 
collected in Table I. 



A. Renormalization of the QEP 

To begin with, we need to regularize the one-loop 
fermion contribution in Eq. ((5J) namely 



T^fcrmion 
^1 — loop 



\QftB\ 
2ir 



IV. RENORMALIZED QM MODEL 

In this Section, we make semi-analytic estimates of the 
polarization and the magnetic susceptibility of the quark 
condensate, as well as for the chiral condensate in mag- 
netic background, within the renormalized QM model. 
This is done with the scope to compare the predictions of 
the renormalized model with those of the effective mod- 
els, in which an ultraviolet cutoff is introduced to mimic 
the QCD effective quark mass. 

In the renormalized model, we allow the effective quark 
mass to be a constant in the whole range of momenta, 
which is different from what happens in QCD (45|. Thus, 
the higher Landau levels give a finite contribution to the 
vacuum chiral condensate even at very strong fields. This 
is easy to understand: the ultraviolet cutoff, A, in the 
renormalized model can be taken larger than any other 
mass scale, in particular A 3> leSl 1 / 2 ; as a consequence, 
the condition p 2 + 2n\eB\ < A 2 is satisfied taking into 
account many Landau levels even at very large eB. The 
contribution of the higher Landau levels, once renormal- 
ized, appears in the physical quantities to which we are 
interested here, in particular in the chiral condensate. 

Since the computation is a little bit lengthy, it is useful 
to anticipate its several steps: firstly we perform regular- 
ization, and then renormalization, of the QEP at zero 
magnetic field (the corrections due to the magnetic field 
turn out to be free of ultraviolet divergences). Secondly, 
we solve analytically the gap equation for the a conden- 
sate in the limit of weak fields, and semi-analytically 
in the opposite limit. The field- induced corrections to 
the QEP and to the solution of t he g ap equation arc 
divergence- free in agreement with [271 ] . and are there- 
fore independent on the renormalization scheme adopted. 
Then, we compute the renormalized and self-consistent 
values of the chiral condensate and of the magnetic mo- 
ment, as a function of eB, using the results for the 
gap equation. Within this theoretical framework, it is 
much more convenient to compute (//) and £/ by tak- 
ing derivatives of the renormalized potential; in fact, the 
computation of the traces of the propagator in the renor- 
malized model is much more involved if compared to the 
situation of the non-renormalized models, since in the for- 
mer a non-perturbative (and non-trivial) renormalization 
procedure of composite local operators is required [54j]. 
Finally, we estimate \i as weu as the behavior of the 
polarization as a function of eB. 



^ f-\-00 77 

/ ^(k* + 2n\Q f eB\+ m ^ 1/2 



71=0 



2tt 



(31) 



To this end, we define the function, V(s), of a complex 
variable, s, as 



V(s) = -N C Y, 



f 



lg/ggj 
2tt 



EM —(k 2 + 2n\Q f eB\+m 2 q y 



(32) 



The function V(s) can be analytically continued to 
s = 0. We define then lf^° n = Hm s ^ + V(s). Af- 
ter elementary integration over k, summation over n and 
taking the limit s — > + , we obtain the result 

Vl~ = N C J2^P^ (- log(2|Q/eB|) + a) B 2 (q) 



-N r 



* — ' 2ir z 



^ \Q f eB\m 2 n [2 
-^E oJ " f--logK) + a) , (33) 



8tt 2 



where we have subtracted terms which do not depend ex- 
plicitly on the condensate. In the above equation, Q (t, q) 
is the Hurwitz zeta function; for Rc(i) > 1 and Re(g) > 0, 
it is defined by the series C(*;ff) = 2~2^Lo( n + ?) _t ; the 
series can be analytically continued to a meromorphic 
function defined in the complex plane t ^ 1 . Moreover we 
have defined q = (to 2 + 2\QjeB\)/2\QfeB\; furthermore, 
a = 1 — 7b— ip(— 1/2), where 7_e is the Eulero-Mascheroni 
number and ip is the digamma function. The derivative 
(' ( — l,q) = d((t,q)/dt is understood to be computed at 
t = -1. 

The first two addenda in Eq. (|3"3"|) arise from the higher 
Landau levels; on the other hand, the last addendum is 
the contribution of the LLL. The function B 2 is the sec- 
ond Bernoulli polynomial; using its explicit form, it is 
easy to show that the divergence in the LLL term in 
Eq. (|33[) is canceled by the analogous divergence in the 
first addendum of the same equation. It is interesting 
that the LLL contribution, which is in principle diver- 
gent, combines with a part of the contribution of the 
higher Landau levels, leading to a finite result. This can 
be interpreted as a renormalization of the LLL contri- 
bution. On the other hand, the remaining part arising 
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Method 


X (GeV- 2 ) 


Ren. Point (GeV) 


Ref. 


Sum rules 


-8.6 ±0.24 


1 


[31] 


Sum rules 


-5.7 


0.5 


[32] 


Sum rules 


-4.4 ±0.4 


1 


[33] 


Sum rules 


-3.15 ±0.3 


1 


[34] 


Sum rules 


-2.85 ±0.5 


1 


[35] 


OPE + Pion Dominance 


-N c /(Av 2 Fl) 


0.5 


[36] 


Holography 


-1.0757V c /(47r 2 F 2 ) 


< 1 


[37] 


Holography 


-iV c /(47r 2 F 2 ) 


< 1 


f38] 


Instanton vacuum 


-2.5 ±0.15 


1 


f39] 


Zero mode of Dirac Operator 


-3.52 


1 


[40] 


Lattice 


-1.547(3) 


2 


[41] 


NJL model 


-4.3 


0.63 


This work 


QM model 


-5.25 


0.56 


This work 



TABLE I. Magnetic susceptibility of the quark condensate obtained within several theoretical approaches. In the table, 
— 130.7 MeV. See the text for more details. 



from the higher Landau levels is still divergent; this di- 
vergence survives in the B — > limit, and is due to the 
usual divergence of the vacuum contribution. We then 
have 



Til { 2 

V£\Z7 = JVcE^ r log(2|Q/eB|) ± a 



I 



\QfeB\rr 



8tt 2 



2 2 

q . rn 2 
log 



2\Q f eB\ 



-N r 



E^c-(-i, 9 ) 



(34) 



The divergence of the vacuum energy is made explicit 
in the expression in the square brackets in Eq. (|3"4")l . The 
mathematical structure of the divergence, namely a pole 
in s = and a logarithm with a dimensional argument, 
is similar to that obtained within the dimensional regu- 
larization scheme. The scale of the logarithm is hidden 
in the 1/s term, and appears explicitly when the diver- 
gence is subtracted. Such a divergence affects only the 
B = effective potential; the corrections due to the mag- 
netic field are either finite or independent on the conden- 
sate. As a matter of fact, in the zero magnetic field limit 
q — >• oo, the fermion bubble becomes 



N r N f 



16tt 2 



log m„ 



1 

a+ 2 



where we have used the relation 



log(g) 



V , 



(35) 



(36) 



with £?2 corresponding to the second Bernoulli polyno- 
mial. Using Eqs. ([3"5[) and ([Ml) we notice that the pole 



2/s is cancelled in the difference V x = V^'oop 1 - V 



Vi 



2\Q f eB\ 




To remove the divergence of the vacuum term, we 
follow Ref. [I?} and adopt the renormalization condi- 
tions that at zero magnetic field, the quantum correc- 
tions do not shift the classical expectation value of the 
a field in the vacuum, (a) = f^, as well as the clas- 
sical value m 2 . To fulfill these conditions, two coun- 
terterms have to be added to the effective potential, 
V c t - — S\ x er 4 /4 ± Sv x <j 2 /2. This amounts to require 



diVo + V^-) 



da 



d 2 (V Q + V c - t -) 



da 2 



= 



(38) 



Taking into account the one-loop divergent contribu- 
tion and the condition (|38p. we can write the renormal- 
ized potential at zero field as 

1 2 

(39) 



V = - (a - v 2 ) - ha : - J — q lop 



16tt 2 



9 2 !l 



where A and v are the renormalized parameters, whose 
expression is not needed here. Equation (|3"9"]) is in agree- 
ment with the textbook result for the one-loop effective 
potential of a linear sigma model coupled to fermions by 
means of a Yuwaka-type interaction [521 ]. This is a con- 
sistency check of our regularization and renormalization 
procedures. 



B. Approximate solutions of the gap equation 

Weak fields. We now restore the magnetic field, and 
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firstly take the weak field limit. In this context, small 
field means eB <C m 2 . Using again Eq. (|3l))) we notice 
that the derivative of the (— function cancels the other 
addenda in Eq. Q37p. and the remaining contribution is 



/ 



= -N B 



f 



24?r 2 
24tt 2 



2\Q f eB\ 



(40) 



which is in agreement with the result of [27| ■ In the above 
equation we have followed the notation of [43[ introduc- 
ing an infrared scale [i, isolating and then subtracting 
the term which does not depend on the condensate. The 
scale /i is arbitrary, and we cannot determine it from 
first principles; on the other hand, it is irrelevant for the 
determination of the a— condensate. We expect fj, « 
since this is the typical scale of chiral symmetry breaking 
in the model for the a field. The correction (|4U)l lowers 
the effective potential of the broken phase, thus favor- 
ing the spontaneous breaking of chiral symmetry. The 
result in Eq. (|4T)|) is UV divergence-free as anticipated; 
thus, it is independent on the particular procedure used 
to regularize the QEP. 

In this limit, it is easy to obtain analytically the behav- 
ior of the constituent quark mass as a function of eB. As 
a matter of fact, we can expand the derivative of the QEP 
with respect to cr, around the solution at B = 0, writing 
(a) = f-rr + 5a. Then, a straightforward evaluation leads 
to 



m q = g.U 1 + - 



9 127T 2 /2 m ; 



-(eB)' 



(41) 



As anticipated, the scale jj, is absent in the solution of 
the gap equation. 

Strong fields. In the limit eB ^> ?« 2 , we can find an 
asymptotic representation of V\ by using the expansion 
£'(— 1, q) = Co + Ci(q— 1) valid for q jw 1, with c = —0.17 
and ci = -0.42. Then we find 



E"<Vr 
8# 



\Q f eB\ log 



2\Q f eB\ 



-N r . 



f 



\QfeB\_m 2 q 
2tt 2 



-ci 



(42) 



where we have subtracted condensate-independent terms. 

In the strong field limit it is not easy to find analyt- 
ically an asymptotic representation for the sigma con- 
densate as a function of eB; therefore we solve the gap 
equation numerically, and then fit data with a convenient 
analytic form as follows: 



bleBl 1 / 2 



eft 
\eB\ 



(43) 



where b = 0.32 and c = 32.78. At large fields the quark 
mass grows as le-Bl 1 ' 2 as expected by dimensional anal- 
ysis; this is a check of the equations that we use. 



C. Evaluation of chiral condensate and magnetic 
moment 

Chiral condensate. To compute the chiral condensate 
we follow a standard procedure: we introduce source 
term for //, namely a bare quark mass mj, then take 
derivative of the effective potential with respect to m/ 
evaluated at rrif = 0. This amounts to derive only 
V^-Yoop 11 ' since both the classic potential and countert- 
erms do not have a dependence on rrif. According to 
Eq. ([37]) we find 



iff) 



dV 



drrif 



dV x 



dm 



(44) 



The first and second addenda on the r.h.s. of Eq. (j44j) 
represent the vacuum and the field-induced contributions 
to the chiral condensate respectively. The vacuum contri- 
bution is divergent, and can be renormalizcd according 
to the procedure of renormalization of composite local 
operators [54j . It is not necessary to perform this pro- 
cedure here, because the numerical value of the vacuum 
condensate is not necessary in our discussiorQ. There- 
fore, it is enough to compute only the contribution at 
B/0 arising from V\, which is finite. 

In particular, for the weak field case we obtain 



(ff) = (/7)o 



N c \Q f eB\ 
12tt 2 m n 



(45) 



On the other hand, in the strong field limit we have 



(ff) 



N c m 
"^ 2 



1 (\QfeB\ 



log 



2\Q f eB\ 



(46) 



Using Equations (|41j) and (j43j) . we show that the chiral 
condensate scales as a + b(eB) 2 for small fields, and as 
| e £|3/2 f or i ar g e fields. 

Magnetic moment. Next we turn to the computation 
of the expectation value of the magnetic moment. The 
expression in terms of Landau levels is given by Eq. (|17D , 
which clearly shows that this quantity has a log-type di- 
vergence. In order to avoid a complicated renormaliza- 
tion procedure of a local composite operator, we notice 
that it is enough to take the minus derivative of V\ with 
respect to B to get magnetization, M. [43| , then multiply 
by 2ml Qf to get the magnetic moment. This procedure 
is very cheap, since the £?— dependent contributions to 
the effective potential are finite, and the resulting expec- 
tation value will turn out to be finite as well (that is, 
already renormalized). 



1 In principle, in the renormalization program of the composite 
operator (qq) one can impose, as a renormalization condition, 
that the chiral condensate in the vacuum is consistent with the 
results of |49| at the renormalization point M = 1 GeV. 
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In the case of weak fields, from Eq. (|4TJ)) we find 



Q f \eB\m q rr± 



6tt 2 



(47) 



On the other hand, in the strong field limit we get from 
Eq. (JUJ) 



N, 



m q 2\Q f eB\ 
log — — 



4tt 2 



(48) 



The above result is in parametric agreement with the 
estimate of magnetization in J43|. In fact, m q ps \eB\ 
in the strong field limit, which leads to a magnetization 
M ps B\ogB. 

Using the expansions for the sigma condensate at small 
and large values of the magnetic field strength, wc argue 



that E/ ps |eS| 
strong field. 



in a weak field, and E / 



\eB\ 3 / 2 



m a 



D. Computation of chiral magnetization and 
polarization 

We can now estimate the magnetic susceptibility of the 
quark condensate and the polarization as a function of 
eB. For the former, we need to know the behavior of 
the magnetic moment for weak fields. From Eq. (|47|) and 
from the definition (TTJ) we read 

The presence of the infrared scale fj. makes the numerical 
estimate of \ uncertain; however, taking for it a value [i ps 
/tt, which is the typical scale of chiral symmetry breaking, 
we have x(ff) ~ 44 McV, which is in agreement with 
the expected value, see Eq. (0. In Fig. [4] we plot /(/x) 
as a function of /i. The interval on the /i-axis delimited 
by the green and the blue vertical lines is the range in 
which we obtain a value of x which is consistent with 
phenomenology 

Next we turn to the polarization. For weak fields we 
find trivially a linear dependence of (if on \QfeB\, with 
slope given by the absolute value of x in Eq. (|49[) . On the 
other hand, in the strong field limit wc find, according to 
Eq. - 



m\ + \Q f eB\ 



1 - 



\Qf\ 



(50) 



where we have used Eq. (|4"o]) . This result shows that 
the polarization saturates at large values of eB, but the 
asymptotic value depends on the flavor charge. 

It is interesting to compare the result of the renormal- 
ized model with that of the effective models considered 
in the previous Section. In the former, the asymptotic 
value of (if is flavor-dependent; in the latter, (if — > 1 
independently on the value of the electric charge. Our 




0.4 0.6 0.8 1 

Mr 



1.2 1.4 



FIG. 4. Magnetic susceptibility of the quark condensate mul- 
tiplied by the chiral condensate, in MeV, as a function of the 
infrared scale (i, as given by Eq. ((49}. The interval on the 
/i-axis delimited by the green and the blue vertical lines is 
the range in which we obtain a value of x which is consistent 
with phenomenology. 



interpretation of this difference is as follows: comparing 
Eq. (|50p with the general model expectation, Eq. , we 
recognize in the factor \Qf\/(b+ \Qf\) the contribution 
of the higher Landau levels at zero temperature, which 
turns out to be finite and non-zero after the renormal- 
ization procedure. This contribution is then transmit- 
ted to the physical quantities that we have computed. 
The trace of the higher Landau levels is implicit in the 
solution of the gap equation in the strong field limit, 
namely the factor b in Eq. (|43p , and explicit in the addi- 
tional \Qf\ dependence in Eq. ([50)) . A posteriori, this 
conclusion seems quite natural, because in the renor- 
malization procedure we assume that the effective quark 
mass is independent on quark momentum, thus there is 
no cut of the large momenta in the gap equation (and 
in the equation for polarization as well). In the effec- 
tive models considered in the first part of this article, 
on the other hand, the cutoff procedure is equivalent 
to have a momentum-dependent effective quark mass, 
m q = ga<d(A 2 — p 2 — 2n\Q feB\), which naturally cuts 
off higher Landau levels when eB ^> A 2 . At the end 
of the days, the expulsion of the higher Landau levels 
from the chiral condensate makes (if 1 in the strong 
field limit. Our expectation is that if we allow the quark 
mass to run with momentum and decay rapidly at large 
momenta, mimicking the effective quark mass of QCD, 
higher Landau levels would be suppressed in the strong 
field limit, and the result ((50)) would tend to the result 
in Fig. d 



V. CONCLUSIONS 

In this article, we have computed the magnetic suscep- 
tibility of the quark condensate, x, and the polarization, 
(if, in a background of a magnetic field, B, by means of 
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two chiral models of QCD: the quark-meson model and 
the Nambu-Jona-Lasinio model. The knowledge of these 
quantities is relevant both theoretically and phenomeno- 
logically. Indeed, the magnetic susceptibility of the quark 
condensate might lead, in the photoproduction of lepton 
pairs, to an interference between the Drell-Yan process 
and the photon-quark coupling, the latter induced by the 
chiral-odd distrubution amplitude of the quark [HJ . The 
two models are widely used to study the phase diagram 
of QCD in several regimes; it has been proved in differ- 
ent contexts, that they offer a good theoretical tool to 
compute low-energy QCD properties. It is thus inter- 
esting to compute, within these models, quantities which 
characterize the QCD vacuum in a magnetic background, 
and compare the results with those obtained within dif- 
ferent theoretical frameworks, both analytically and nu- 
merically. 

In the first part of this article, we have reported our 
results about x an d fJ- / obtained within a numerical self- 
consistent solution of the model, in the one-loop approxi- 
mation. Our results on \ are summarized in Table I, and 
are in fair agreement with previous estimates. Besides, 
our data on polarization are collected in Fig. 3. For the 
latter, we obtain a saturation to the asymptotic value 
/ioo = 1 at large values of eB, which is understood as a 
lowest Landau level dominance in the quark condensate. 
This is in agreement with recent results obtained within 
Lattice QCD simulations with two colors and quenched 
fermions (4f| . 

In the second part of the article, we have estimated \ 
and /if within a renormalized version of the QM model. 
The numerical value of \ m this case is quite uncertain 
because of the presence, in the final result, of an unknown 
infrared scale \i. However, as shown in Fig. 01 taking for 
\i a numerical value around f n , which is the typical scale 
of chiral symmetry breaking in this model, we obtain a 
value of x which is consistent with phenomenology. 

Within the renormalized model we find a saturation of 
fif at large B, in qualitative agreement with our findings 
within the effective models in Section III and quenched 



QCD 4l|. In the renormalized model, the asymptotic 
value of fif is flavor-dependent; in the effective models 
and quenched QCD, (if — > 1 independently on the value 
of the electric charge. We attribute this difference to the 
presence, in the theory, of the renormalized contribution 
of the higher Landau levels: in the renormalized model, 
the higher Landau levels give a finite contribution to the 
chiral condensate at zero temperature, even in the case of 
very strong fields. After renormalization of the QEP, this 
contribution is finite and non-zero, and is transmitted to 
the physical quantities that we have computed. 

The results we obtained are quite encouraging, and 
suggest that a systematic study of external field induced 
condensates at zero, as well as finite, temperature within 
chiral models is worth to be done. As a natural contin- 
uation of this study, it would be interesting to add the 
strange quark, following the work [39[. Moreover, Lattice 
simulations have shown that x is almost insensitive to the 
temperature, at least in the confinement phase (4l| . This 
result is achieved with quenched fermions. On the other 
hand, the chiral model handle with dynamical fermions; 
it would be therefore of interest to address the question 
of the behavior of x as a function of temperature using 
chiral models. Furthermore, it has been shown [53j that 
local fluctuations of topological charge induce a quark 
electric dipole moment along the direction of a strong 
magnetic field. This problem can be studied easily within 
the chiral models, introducing a pseudo-chemical poten- 
tial -Up conjugated to chirality imbalance as already done 
in 0, Ill- 
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